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Abstract. We investigate the transformation of the distribu- 
tion function in the relativistic case, a problem of interest 
in plasma when particles with high (relativistic) velocities 
come into play as for instance in radiation belt physics, in 
the electron-cyclotron maser radiation theory, in the vicin- 
ity of high-Mach number shocks where particles are acceler- 
ated to high speeds, and generally in solar and astrophysical 
plasmas. We show that the phase-space volume element is 
a Lorentz constant and construct the general particle distri- 
bution function from first principles. Application to thermal 
equilibrium lets us derive a modified version of the isotropic 
relativistic thermal distribution, the modified Jiittner distribu- 
tion corrected for the Lorentz-invariant phase-space volume 
element. Finally, we discuss the relativistic modification of a 
number of plasma parameters. 

Keywords. Relativistic invariance of the one-particle distri- 
bution function, Jiittner distribution, relativistic plasma pa- 
rameters 



1 Introduction 

One-particle distribution functions are defined as the proba- 
bility of finding a certain number of particles of species 
a at time f in a given phase-space volume element dx^dp^ as 



/a(x,P,0 



dx^dp^ 



(1) 



Since the above function is a probability function this expres- 
sion can also be read as a normalisation condition 



(2) 



saying that the product of the probability divided by the num- 
ber density of particles in the phase space volume under con- 
sideration must be one, thus being invariant. Under special- 
relativistic conditions the definition of such a distribution 
poses a non-trivial problem. In this Letter we define such a 
distribution from first principles showing that of the critical 
points the first is to confirm that the phase-space volume el- 
ement must correctly be made Lorentz invariant. Once this 
has been done, it is possible to define a distribution function 
valid for the entire phase space. 

2 Construction of phase-space distributions 

Clearly, if the distribution function is a proper probability, 
it is by its nature invariant with respect to any coordinate 
transformation of the spatial and momentum space volume 
elements, if the number of particles in the volume element 
is kept constant. This follows from Liouville's theorem and 
the dynamics of the particles (at least as long as no parti- 
cles are lost or added to the volume element under consid- 
eration). However, a proper transformation of the volume 
element might be in place. In the Galileian non-relativistic 
case there is no problem at all. Problems arise when the 
particle dynamics becomes relativistic. While the functional 
form of the phase-space probability in such transformations 
is conserved, its dependence on the new phase-space coordi- 
nates x',p' will change during the transformation. This has 
been noticed long ago (fo r the basic arguments see, e.g.. 
Landau and Lifshitzl Il975[ Section 10), and the transforma- 
tion of the phase-space volume becomes non-trivial. In fact, 
in order to arrive at a consistent formalism and understanding 
one cannot simply transform the above formula; rather one 
must return to the exact definition of the distribution func- 
tion given by 
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where x,(f),p,(f) are the dynamical coordinates of the /th par- 
ticle. Since F sums up all points in phase space which are 
occupied by the ^ particles numbered /=!,..., it is ex- 
act. Accordingly, the exact distribution of one single particle 
is defined as 

Fj{x,p) = 8[x-xj{t)]8[p-pjit)] (4) 

and the exact distribution of a group 8^ of particles is given 
by 

8F,..j(x,p) = ^^Yl8[x~Xi{t)]8[p~pdt)] (5) 

This all looks simple and is simple as long as one does not 
include relativity. Then, however, particles are all in motion, 
each in its own proper frame and with its own proper di- 
rection of velocity and momentum. The transformation from 
that frame into another frame, i.e. into the observer's station- 
ary frame involves the Lorentz transformation of each parti- 
cle's frame separately. In the (primed) particle's proper frame 
(x',p') the particle is located at the origin and the proper- 
frame distribution f Eq. (HI reads 

f0(x',p') = 5(x')5(p') (6) 

while Eq. (|4]i is the expression in another frame with re- 
spect to which the particle is displaced, has momentum P; (r) 
and moves with relative velocity \j{xj,t) with respect to the 
new frame. Taking the axes of the two frames all parallel 
to each other and assuming that the particle instantaneously 
moves along x', and defining j3=v/c,7= 1/ \/l — = 
a/I + /rrflc^ with m particle mass, momentum p = mcyP, 
thus /3 = 7^ ' -^72 — I, and time t^ct one has as usually for 
the spatial coordinates 

[x', y, z', t'] = [y{p){x~Pt), y, z, 7{p)it-M] (7) 

Inserting into Eq. (|4|i then yields for the transformed exact 
single-particle distribution function 

Fjix,p,t) = 8{x- [Yj{p){x-pjt), y,z]]x 

X 8[p-p'jit)]8[t- 7j{p){t^l5jx)] (8) 

In this expression p^- = pj because there is only one (instan- 
teneous) relative velocity between the moving jth particle's 
and the observer's frame. However, j3 is a function of 7(75), 
and though the functional form of the probability distribu- 
tion (represented by the product of 5-functions) is conserved, 
the dependence on phase-space coordinates becomes compli- 
cated. 

In considering many particles with their different proper 
velocities, proper momenta and different directions of veloc- 
ities, one needs to rotate the axes of all the proper frames 
in performing the various Lorentz transformations. This in- 
hibits the use of the exact equation Eq. (|5]l even in the case 



of a completely isotropic angular distribution of the various 
velocities. Moreover, stationarity of the distribution becomes 
illusionary for several reasons, the simplest one that it is im- 
possible to assign initial conditions to all particles at a fixed 
time f = for the obvious reason that this requires time. This 
time is needed for jumping from one proper frame to the next 
and, thus, is limited by relativity and the requirements of rel- 
ativistic synchronisation. It can be achieved only for a small 
limited number of particles. 

3 Lorentz invariance of full phase-space volume ele- 
ment 

The only way out of this dilemma is to neglect the proper 
particle motions altogether and assume that by some myste- 
rious (thermodynamic) reason the particles have managed to 
reach a stationary final state where their momenta are statisti- 
cally distributed. In this case a common phase-space volume 
can be attributed to them all, and each particle moves with 
its (thermal) momentum with respect to this common proper 
frame. If the momenta are isotropically distributed in space 
then the probability distribution will not depend on the direc- 
tion of velocities. It can then be considered as averaged over 
all angles and depends just on the modulus of the momentum 
P- 

The way to do this is to take one single particle and to 
smooth it out so as to fill a certain infinitesimal but finite 
proper volume d^x'jd^pp which is actually a volume ele- 
ment, while defining a probability of finding the particle 
in this volume fj{x',p') which now is a continuous func- 
tion which vanishes outside the volume element and is nor- 
malised according to 

l=lfj{x'j,p'j)d%d''pj (9) 

with integration just over the infinitesimal 7th volume. It 
is unity because the infinitesimal volume contains only the 
one yth particle, and fj = outside it. In principle the latter 
fact would allow for extending the integration over the entire 
phase space volume. 

Within this small volume element it can be said that pre- 
scribing the initial state has been done at a fixed time t' = 0. 
This is correct for sufficiently small volumes such that syn- 
chronism is warranted. It thus becomes possible to Lorentz 
transform the infinitesimal volume to the un-primed observer 
frame. Since the particle velocity at fixed time is linear, one 
may rotate both frames to coincide in the directions of their 
axes. 

Applying the spatial transformation Eq. dTji at f = and ob- 
serving that only the coordinate x' is affected by the Lorentz 
transformation becoming contracted when transforming to 
the un-primed observer frame, one immediately finds (cf., 
also, .Weinberg. ,1973l his Eqs. (2. 1 0. 1 1 ) & (2. 1 0. 1 2) for the 
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transformation of the particle 4-current density) that 



(10) 



where the index j at d Xj reminds one that the observer's 
frame has been rotated into the particle frame, and py is the 
particle momentum of the smeared out infinitesimal volume 
of the particle in the observer's frame. In the proper frame 
of the particle the momentum has zero average (p^) ~ 0. 
In the frame of the observer the finite infinitesimal momen- 
tum space volume d^pj is located at the tip of the particu- 
lar particle-momentum vector py. The proper Lorentz factor 
Yj{p) depends only on the modulus of py. 

One may also understand the small infinitesimal momen- 
tum space volume as the thermal uncertainty of the particle 
momentum that is caused by the quiver velocity of the parti- 
cle. This quiver velocity is the rms velocity {u)j = \J (Auy)^ 
of the particle's proper motion, its thermal speed. It may be 
assumed to be isotropic though this is not necessary if this 
thermal speed is anisotropic due to some other reason. 

With all this in mind we can proceed to the transformation 
of the momentum space element. This is more involved than 
the transformation of the spatial volume. It requires return- 
ing to the definition of the momentum py = mcYj{p)Pj and 
taking into account that the particle energy and momenta in 
relativity are not independent, being related via the relation 
pc = S" ox S ~ mY[p)c^. This implies that the transforma- 
tion of momentum and energy has to be done together This 
becomes 



" S' ' 


= rip) 


■ 1 -P' 




■ s ' 


cp'. 






CPx 



p'. 



Py 

Pz 



The variations in energy and momentum become 



(11) 

(12) 
(13) 



dS' = -y{p)[c^dpx-dS] =0 

The last expression follows because, in the particle rest 
frame, the energy remains unaffected by the transformation, 
yielding as usually dS" = cfidpx- Hence, 

dpx 



dp',^yip)dpx{l-l5^) = 



7{p) 



(14) 



which, of course, is just a special case of the well-known 
transformation law 



dpx 



dp'x 

inc'Y{p') mcjip) 



, with y{p') = 1 



(15) 



and the single-particle phase-space volume including Eq. 
(flOb is an invariant under Lorentz transformation 



d^x'jd^p'j = d^Xjd^pj 



(16) 



This holds for each particle separately such that one can de- 
fine a single-particle distribution function in the observer's 
frame by 

fj{xj,9j)d\jd^Pj = fj{^j,9))d\'/p'j (17) 



4 Global phase-space distribution 

One would wish to add up all single particle distributions to 
one distribution which covers the entire phase space. This 
task encounters the problem that so far the two frames, the 
y th particle frame and the observer frame, have been aligned. 
In order to obtain a similar result this must be done for all 
particles. Hence, simply adding the distribution functions is 
inappropriate. One needs to introduce a distribution of angles 
for the various particle moments as seen from the stationary 
observer frame. This can be achieved with the help of the 
solid angle £2p and the distribution <J>(np). Then in the sta- 
tionary frame one defines the probability of finding all ^ 
particles in the volume of phase space occupied by all the 
smeared-out infinitesimal though finite volumes of the parti- 
cles (assuming that they do not overlap, i.e. in the absence of 
collisions) 

.A' 

/(x,p) d\ d^pdQ.p = n//(x;,p;)<I'(np) dQ.p d% d'p'j 

- n/;[4(X'P)'P;(P''^)]'J'(^p) <^^P d'p (18) 

7=1 

For an isotropic distribution of particle momenta the integral 
over the solid angle gives just a factor An. This yields for the 
probability distribution in the stationary frame the final result 

/(X,p,np) dD.p = n/7[X;(x,P),P ;(P,'^')]«>(np) dD.p 

7=1 

(19) 

Normalisation of this distribution must as usually be done to 
the total number jV of particles respectively the average par- 
ticle density = j'Y . With this definition, the relativistic 
particle phase-space distribution can indeed be written as a 
single average distribution while the directions of the differ- 
ent single particle velocities and thus the various rotations 
of the JV proper coordinate frames is taken care in the an- 
gular dependence of /(x,p,np), which eliminates the oth- 
erwise required dependence of the volume element on the 



proper particle-speed direction as given by iDebbasch et al 
(l200l[ and references therein). 



Clearly, this angular distribution must be introduced by 
hand, making a reasonable guess. This, however does not 
provide any problem as it can be chosen to be inherent to 
the initial model. Again, in the case of an isotropic angular 
distribution it just gives, after integration, a factor of \n. 

5 Applications 

The evolution of the relativistic distribution function 
/(XjPjHp) in a dynamical system is subject to the relativistic 
Vlasov equation 



dt 



- [[jr,/]] = 0, (with the Hamiltonian) (20) 
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written in the 3+1-spUt non-covariant representation, where 
[[. . .]] is the classical Poisson bracket. In fact, here ^(x,p) is 
the 1 -particle Hamiltonian. Under equilibrium conditions the 
time dependence is dropped, and in thermal equilibrium the 
left-hand side is identically zero. 

Relativistic tliermal equilibrium distribution 

The stationary thermal equilibrium distribution assumes the 
form of a relativistic variant of the Maxwell-Boltzmann dis- 
tribu tion. Its determinat ion poses a nontrivial problem (cf., 
e.g., Dunkel etal 2007 . and references therein) as the ther- 
mal equilibrium under relativistic conditions is defined more 
subtly, yielding different forms of the so-called Maxwell- 
Jiittner distribution. Its functional form must in fact be de- 
rived from putting the Boltzmann collision term 

0= f d^p'a,ouda„,n [(/'/)*-/'/] x 



X ^{v'-9?+m^c^[y'^{p')Y''{p)~\]/mcY[p') (21) 

to zero which, in the above equation, has been taken to van- 
ish identically. Here C7£„// is the collisional cross section, 
and dD.coii the solid collisional deflection angle, and primed 
and un-primed quantities refer to the two colliding particles, 
while the asterisk refers to quantities after collision. This then 
yields [(/'/)* ~ f f] ~^ ™d, ultimately, the famous expo- 
nential Boltzmann-Gibbs dependence on energy; 



/ = Aexp 



-@Qmc^y{p)-Y^@icp' 



(22) 



The Lagrangean multipliers A,0o,0, must be determined as 
usual from particle number and global energy-momentum 
conservation; the 4-vector components ©y have dimension 
of inverse energy which suggests that they relate to a partic- 
ular inverse-temperature 4-vector Consistency with the col- 
lisionless part of the Vlasov equation yields that A = const 
does not depend on space, while 0v could become a linear 
4-vector function of space. Under isotropic conditions, how- 
ever, only the average energy comes into play, and one may 
put all ©; = 0, / = 1 , 2, 3, keeping only the constant ©o ^ 0. 

General form of distribution. This approach then leads to the 
canonical form of the isotropic relativistic distribution func- 
tion 

frei[7{p)] = A exp [^©Qmc^Yip)] (23) 

with the only remaining problem to determine the two La- 
grangean multipliers A and ©q. The latter can simply be 
guessed from the requirement that the argument of the ex- 
ponential should be dimen sionless, whic h yields ©o = \/T 



(for rigorous argument, cf. Ilsraell 119631 ) identifying it with 



the inverse of temperature T of the macroscopic system of 
particles. Determination of A requires solving for the nor- 
malisation of the distribution. Here the different approaches 
diverge, and there has no unambiguous consensus be found. 



It is argued ( Dunkel et all 2007 ) that, if rigorous Lorentz 
invariance is imposed (as has been used in our strict 
definition of the general non-stationary distribution func- 
tior{3) the correct (non-angular-dependent part of the) rela- 
tivistic thermal-equilibrium distribution should become the 
modified-Jiittnei distribution. (The ordinary Maxwell-Jiittner 



12J 

distribution function was derived by F. Ijiittnerl 1191 ll who 



obtained it imposing translational invariance in momentum 
sp ace only.) 



Dunkel et a]l (l2007h chose a general approach starting from 



the maximum-entropy principle for the Boltzmann-Shannon 
(relative) entropy density (actually the entropy density per 
particle) sq as a functional of the distribution function 



^o[/(pK(p)] = - / A /(P)log [(mc)''/(p) 



(24) 



where i{p) ~ [mc)^ is the constant momentum space den- 
sity needed for normalising the argument of the logarithm, 
and Yp the momentum space volume. This is subject to the 
constraints 



1 



d'pfip), 



i 



d'f{p)e{p) (25) 



where is the average energy density, i.e. the scalar pres- 
sure which is the trace of the energy-momentum tensor T^^ , 
and e(p) ~mc^7(p) the relativistic particle energy. Defin- 
ing the Haar measure to momentum space, they propose the 
modified-Jiittner distribution for arbitrary dimensions d. Un- 
fortunately, however, Lorentz invariance is applied only to 
momentum space. 

On using the Boltzma nn collision-integral approach, 
Chacon- Acosta et al (l2O10l and references therein) obtain 
another form in arbitrary dimensions d which, for c/ = 3, re- 
duces to t he original M axwell-Jiittner function (already given 
earlier by Ilsraell 119631) . 



JUttner-like thermal distribution. However, all these ap- 
proaches only consider the momentum-space part of the vol- 
ume element neglecting the transformation of the spatial co- 
ordinates, thus keeping the factor = mcy[p) in the denom- 
inator of the momentum-space volume element. This is either 
not justified at all or it is argued that the particles are all con- 
fined to a fixed box which is unaffected by the Lorentz trans- 
formation and invariance. However, the momentum and con- 
figuration space volume elements the product of which forms 
the phase-space volume element, are not independent, as we 
have demonstrated above. Even in this case of a fixed outer 
box, the particle's proper spaces experience linear Lorentz 
contractions when seen from the stationary frame of the ob- 
server, i.e. from the box-frame perspective. The consequence 



'it is important to note that this is the only case which accounts 
for the required general Lorentz invariance of the complete phase- 
space volume element, and only under the condition of its validity 
is it possible to define a global phase-space distribution function. 
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is that the extra proper Lorentz factor jip) in the phase-space 
volume element cancels thereby guaranteeing and restoring 
Lorentz invariance. The correct form of the normalisation 
condition is thus 

A^^= / d'pp'Uiirip)], T^"^ I d'ppWeiirip)] (26) 



where A^^ is the particle-current density 4- vector, and Tf^^ 
the average energy-momentum tensor When taking into ac- 
count the functional dependence of the relativistic distribu- 
tion Eq. ( |23] l on 0, both can be derived from the Integral 



d'pfreliYip)] 



(27) 



by differentiation with respect to 0. For instance, with = 
mcy, one has for the average particle density in the observer 
frame (the stationary box) 



n=n" = -- 



(28) 



Calculatio n of the integral is straightforward. Following ear- 
lier w ork (Juttner 1911 : Israeli 1963 Chacon- Acosta et al 
one introduces spherical-harmonic coordinates in mo- 
mentum space and uses the usual representation of the 
Lorentz factor by hyperbolic functions as 



p^ = OTC cosh^ 

p' = OTC sinh^ (sin0cos(/),sin0sin0,cos0) 
with / = 1,2,3 and < ^ < oo. This yields 

= 47r(OTc)3 / t/^sinh^^cosh^e-®'""''''^"^''^ 



(29) 



(30) 



If one choses x = cosh^ as dummy integration variable, the 
boundaries of the integral become 1 <x<°o. Performing the 
integration yields 



J' = An{mc 



P^^(ot£20o) 

OTc20O 



00 = ^ 



(31) 



an expression in which the second order modified Bessel 
function ^'2(2) appears instead of K\{z). This result is dif- 
ferent from those obtained without taking into account the 
correct Lorentz invariance of the full phase-space volume el- 
ement, including the configuration space. 

Taking the derivative as prescribed above we find 



A = 



2 [ I y-T 1 I y-T 



(32) 



This is slightly more complicated than Jiittner's normalisa- 
tion but takes into account the correct Lorentz invariance of 
the phase-space volume element as required by relativity. In 
any application it is therefore preferable over the former. Eqs. 



( I23I 1, ( |32] |. and ( |3T1 ) define the correct relativistic isotropic 
thermal equilibrium distribution. 

Like the conventional Juttner distribution, the modified 
Juttner distribution given here refers to thermal equilibria 
being just the relativistic version of the Boltzmann energy 
distribution when taking into account Lorentz invariance of 
the phase-space volume element. Observations of high en- 
ergy particles, for instance Cosmic Ray spectra or spectra 
of energetic particles in the solar wind and in the vicinity 
of shocks, frequently deviate quite strongly from such distri- 
butions exhibiting long non-thermal power law tails in en- 
ergy and velocity. The physics of generation of such tails 
is still a problem of intense research. Clearly such distri- 
butions are not described by thermal equilibria; it is usually 
assumed that highly non-linear wave particle interactions in 
non-equilibrium systems, i.e. mostly in the presence of col- 
lisionless shock waves, or else plasma turbulence, would be 
capable of generating particle spectra and distributions with 
extended tails. Since in most cases those tails reach far into 
the relativistic domain one would from the beginning as- 
sume that the non-linear problem of tail generation should 
be treated relativistically working with relativistic thermal 
distributions as initial equilibrium conditions which, when 
exposed to turbulence or non-linear processes, finally evolve 
relativistic tails. These questions are particularly important 
in astrophysical applications to energetic particle spectra. 

There have also been attempts to describe power law 
tail distribution, called if-distributions, as some intermediate 
thermodnamic equilibrium in non-linear or turbulent interac- 
tion. Such attempts have so far been fairly artificial and aca- 
demic being either based on a modification of Boltzmann's 
collision integral or on the proposition of some special form 
for the thermodynamic entropy depending on some free pa- 
rameter like K. Distributions of this kind fit fairly well ob- 
servations. In how far they obey a thermodynamic state has, 
however, not been clarified from basic principles. It has also 
not been shown, how the free parameter K can be related to 
internal non-linear plasma processes on which it must neces- 
sarily depend if being real. In principle, of course, one could 
extend formally the derivation of an equivalent to the Juttner 
function also for such distributions. However, as long as there 
is no thermodynamic basis for them this remains a purely 
academic problem. 

Various plasma quantities 

The Lorentz-invariant property of the phase-space volume 
element has consequences for the definition of a number 
of phase-space averages of plasma parameters. These come 
into play via the relativistic variation of the particle mass 
m — > my{p) which depends on the single particle momentum. 
Hence, in presence of large particle numbers the collective 
averages of some quantities imply the phase-space integra- 
tion over all contributing particles. 
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Plasma frequency 

The square of the plasma frequency of particles of rest mass 
rriQa and charge qa is defined as O)^ = q\iNal £Qnia, where 
wJa = fndajip) and is the particle density which follows 
from the zeroth moment of the relativistic distribution func- 
tion. Under relativistic conditions with the single particles 
having proper Lorentz factors 7(/?), this becomes 



d^xd^p 2 
/a(p)^< 



£omoaJ r(p) J Y(p) 

The extra Lorentz factor in the denominator of the integrand 
arises from the relativistic mass-dependence on the proper 
yip). When introducing the rest-frame density A^^ and as- 
suming that the relativistic distribution is homogeneous in 
configuration space the integration over space can be per- 
formed leading to the second form in the above equation, 
where = qa^Oa/^Q'noa is the square of the nonrelativis- 
tic plasma frequency. The above expression takes into ac- 
count that each particle has its own relativistic mass which 
depends on its proper momentum. 

So far we assumed that the plasma in the observer's frame 
is at rest, in which case only the proper motions of the par- 
ticles count. If there is a bulk flow with bulk momentum po, 
then the un-primed frame can be considered as the rest frame 
of the plasma as a whole, and a new observer frame can be 
defined in which the plasma moves with bulk Lorentz factor 
r{po). In this case the average density and mass transform 
like N ^ Nr{pQ) and {ina) (ma)r(po), where from Eq. 

d^p 



(ma) =moa 



(34) 



yielding that the ratio of NT{pq) / {ma)T{pQ) remains inde- 
pendent of r(/7o). Thus, for bulk relativistic flows the bulk 
Lorentz factor T{pq) drops out of the plasma frequency, mak- 
ing it independent of any relativistic bulk flows. 

Relativistic length scales 

The above definition of the plasma frequency leads to the 
relativistic redefinition of some lengths which depend on the 
plasma frequency. These are the Debye length and the inertial 
scales. The Debye length Xd does not depend on the particle 
mass. Because of this reason it remains invariant as long as 
no bulk flow is considered. In a medium flowing with bulk 
momentum po having a bulk Lorentz factor r(/:'o) its density 
dependence causes a weak dependence 

Ad^Aoo/VTX^ (35) 

As for the inertial scales = c/Oa, on the other hand, it 
is clear that the dependence of the plasma frequency makes 
them sensitive to the internal Lorentz factor as contained in 
the definition Eq. (|33l l of the plasma frequency. It is obtained 
by averaging the inverse inertial length taking into account 
that fa is normalised 



(A, 



(36) 



which retains the functional form of the inertial length even 
in the relativistic case yielding for the root.mean-square of 
the inertial length of ath particle species A^"' ~ \J (A^) = 
c/-\/ (ft)^), simply replacing the non-relativistic plasma fre- 
quency with its relativistic average. This form remains un- 
changed even in the case when the plasma moves at bulk 
Lorentz factor r(/:)o), for the reasons explained above. 

Alfven speed, Mach number, magnetisation parameter 
As a last quantity we consider the relativistic Alfven velocity. 
The Alfven speed is defined ?&Va^B/ ^ptomN showing that 
it depends on the average particle mass density mN. Since 
is given as the zero moment of the distribution function, one 
generalises the mass density as 



(37) 



maNa) = / d p y{p)fa{v) 



One may note that particles of high internal momentum con- 
tribute substantially to the relativistic mass density by be- 
coming relativistically heavier Moreover, if the plasma ex- 
hibits a bulk flow with Lorentz factor r(/:'o), then the mass 
density increases as T^. 

Calculation of the relativistic Alfven velocity can be done 
in analogy to calculating the inertial scale by taking and 
averaging its inverse over the phase space. This yields 



B 



^/lM^{maNa) 



(38) 



Here the relativistic particles reduce the Alfven speed due 
to their increasing mass effect. This also implies that the 
Mach number . # = V /Va will increase, if the plasma flows 
with bulk velocity V = cF^ ' y/F^ — 1, in which case the Mach 
number becomes proportional to sa vT^ — 1, where 
^0 is the non-relativistic Mach number 

Finally, a comment on the magnetisation parameter (in- 
verse plasma-j3) is in place. This is defined as the ratio of the 
magnetic to kinetic energy densities 



2pk)NK 



(39) 



where K=T -\- Sj^iu is the sum of thermal and kinetic plasma 
energies, the latter playing a role only in plasma in bulk mo- 
tion. A plasma not in bulk motion has thus constant magneti- 
sation parameter a„i = const. However, when the plasma is in 
bulk relativistic motion, one has N ^ NF{po), and the mag- 
netisation parameter decreases. Moreover, if the bulk energy 
is substantially larger than the thermal energy, as is the case 
in ultra-relativistic cool flows which are important in astro- 
physics, one has oc F{po) ^ T and thus 



B' 



1 



2^NK r2(po) 



<1 



(40) 



which usually is very small and in this form is commonly 
used in ultra-relativistic astrophysical flows. 
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Under this condition a flow can be considered as very 
weakly magnetised unless it becomes capable of generating 
its own strong equi-partition magne tic fie l d. This is th e case 
when the Weibel instability Ci Weibell [l959t iFriedi [l959l) takes 
over under some well defined conditions. 



6 Conclusions 

In summary, we have given here an argument for the Lorentz 
invariance of the phase-space volume element and con- 
structed the corresponding relativistic distribution function. 
This has consequences for the calculation of some relevant 
plasma parameters in the relativistic domain. Application of 
such parameters can be found in the cyclotron-maser theory 
as for instance in the aurora or solar physics as well as in as- 
trophysics, and to radiation belt problems where relativistic 
particles come into play. 

We have also shown that the isotropic thermal equilibrium 
distribution function in a relativistic plasma Eq. (|23] l can be 
given analytically [see Eq. [32l; it has the form of a modi- 
fied Jiittner distribution. Calculation of the anisotropic ther- 
mal distributio n which has been argued is more correct than 
the is othermal ( Israel , 1963 ; van KampenL 1968 ; Nakamura! 
2009h in the relativistic case would be more involved and 
becomes possible only when the 4-vector of inverse temper- 
ature is properly taken into account. Since these questions 
have not yet been clarified even in basic relativistic thermo- 
dynamics, we have not invested any effort into the determi- 
nation of the general relativistic distribution in that case. The 
normalisation condition changes when the exact Lorentz in- 
variance of the phase space volume is taken into account. 
This has consequences for the relativistic ideal gas laws, 
which have not been explored in this letter. The interesting 
temperature-4-vector as well as anisotropic cases pose fur- 
ther hurdles and have not been considered here. 
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